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Regenerative cooling with fuel as the coolant is used in the scramjet engine. In order to grasp the dynamic cha- 


racteristics of engine fuel supply processes, this article studies the dynamic characteristics of hydrocarbon fuel 


within the channel. A one-dimensional dynamic model was proved, the thermal energy storage effect, fuel volume 


effect and chemical dynamic effect have been considered in the model, the ordinary differential equations were 


solved using a 4th order Runge-Kutta method. The precision of the model was validated by three groups of expe- 


rimental data. The effects of input signal, working condition, tube size on the dynamic characteristics of pressure, 


flow rate, temperature have been simulated. It is found that cracking reaction increased the compressibility of the 


fuel pyrolysis mixture and lead to longer responding time of outlet flow. The responding time of outlet flow can 


reach 3s when tube is 5m long which will greatly influence the control performance of the engine thrust system. 


Meanwhile, when the inlet flow rate appears the step change, the inlet pressure leads to overshoot, the overshoot 


can reach as much as 100%, such highly transient impulse will result in detrimental effect on fuel pump. 
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Introduction 


Hydrocarbon-fueled scramjet has broad application 
prospect in the field of space transportation and military 
in the foreseeable future!!! The hydrocarbon carried by 
scramjet not only works as fuel for combustion, but also 
works as coolant'***!, As hydrocarbon fuel flows through 
cooling channels of regenerative cooling system for 
scramjet, there is a dramatic variation in properties of 
hydrocarbon fuel, since both fuel temperature and pres- 
sure would change substantially, especially under the 
operating conditions with endothermic chemical reaction 
of hydrocarbon fuel!®”!. 

As computational fluid dynamics (CFD) is becoming a 
powerful simulation tool'*!, there has been significant 
interest in developing mathematical models for cooling 
channel, in an order to analyze the complicated interac- 
tions between the various factors and to optimize the 
system performance. Most of them focused on the steady 
characteristics!°''!, T.A. Ward studied the pressure ef- 
fects on flowing mildly cracked n-Decane at steady con- 
dition”. Silong Zhang introduced 3-D model of cracking 
coolant in cooling channel, and studied the thermal beha- 


vior in the cracking reaction zone at different aspect ra- 
tios”, Ward et al. conducted studies about the pressure 
effects on flowing mildly-cracked n-Decane using two 
dimensional modell", 

Other researchers investigated the transient dynamic 
characteristics!'*!*!*'°!, Nicolas Gascoin developed a 
one-dimensional transient numerical model using a com- 
plete validated pyrolysis mechanism for n-dodecane!'?!, 
Dahm K. D et al uses a full n-dodecane pyrolysis me- 
chanism (1185 reactions and 153 species, designed by the 
French laboratory DCPR)'"*!, However, each calculation 
costs too much hours and computing resource. This high- 
lights the fact that CFD cannot be suitable for complete 
transient simulation with chemistry effects, particularly 
as an engineering tool on a single computer. Wen Bao 
conducted a one-dimensional unsteady model of high- 
temperature fuel supply tube containing chemical reac- 
tion of hydrcarobon fuel, and the model can reflect not 
only the dramatic variation in properties of hydrocarbon 
fuel caused by the changes of temperature and pressure in 
wide ranges, but also the coupling among flow, heat 
transfer and chemical reaction of hydrocarbon fuel’! Its 
simplified chemical mechanism greatly save the time to 


run a calculation without losing the capacity of analyze 
the flow and thermal characteristics of the cooling chan- 
nel. Although various kinds of model has been built, 
none of the researchers has analyzed the transient dy- 
namic characteristics of hydrocarbon fuel within the 
channel, and their effects on the fuel supply system and 
thermal protection has not been talked about. 

The purpose of this paper is to study the dynamic cha- 
racteristics of cooling channel during the engine working 
process, the one-dimensional unsteady model of high- 
temperature fuel supply tube containing chemical reac- 
tion of hydrcarobon fuel in the paper [15] is used. The 
effects of input signal, working condition, tube size on 
the dynamic characteristics of pressure, flow rate, tem- 
perature have been simulated. The influence of transient 
dynamic characteristics on the engine fuel supply system 
and thermal protection has been discussed. 


Model and solving method 


The model of hydrocarbon fuel within channel under 
supercritical pressures consists of two parts: one is the 
fluid flow in cooling channel, and the other is the proper- 
ties of hydrocarbon fuel. 


Geometry model of cooling channel 


A kind of cooling channel with rectangular cross sec- 
tion is usually used for hypersonic application, and Fig- 
ure 1 shows the corresponding schematic diagram of a 
cooling channel, into which the heat flux that comes from 
engine is injected. Assuming that it is a 1D flow for coo- 
lant in the cooling channel and all the variations in the 
radial fuel properties and in the radial flow velocities can 
be neglected, by using a hydraulic diameter d to charac- 
terize the cooling channel, the model of cooling channel 
is presented in detail below. 


Fig. 1 Schematic diagram of cooling channels with heat flux 
injected into 


Modeling of fluid flow in cooling channel 


The set of PDE (partial differential equations) used to 
describe the 1D flow dynamics in a cooling channel is 
established by applying the conservations of mass, mo- 
ment, and energy to an infinitesimal control volume. The 
mass conservation equation can be expressed as equation 
(1) 

ee + om =0 (1) 
Ot Ox 
The heat generated by the viscous effect of fluid, axial 


heat conduction of fluid, and kinetic energy of fluid are 
usually much smaller than the radial heat which is in- 
jected into the fluid. Then, by neglecting those relatively 
smaller terms, the energy conservation equation can be 
expressed as equation (2) 
ð fa] 4 
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where qr is the convective heat between the fluid and 
the inner wall surface of the cooling channel. It can be 
expressed as equation (3) 
qs =a-(T,-Tp) (3) 
where œ is a heat transfer coefficient determined by 
equations (4) and (5). Here, only a qualitative expression 
is given, and the specific heat transfer correlation em- 
ployed in this study will be discussed later on. Tw is the 
wall temperature of the cooling channel, and it is calcu- 
lated through equation (6) by neglecting its radial tem- 
perature difference. Although axial heat conduction for 
the wall of the cooling channel is neglected, the axial 
wall temperature is still varying; so, expression O7,/0¢, 
instead of dT,,/dt is used in equation (6) 
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Re and Pr used in equation (5) are expressed as equa- 
tions (7) and (8), where vector u has become a scalar 
through the absolute value transformation 
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Considering that cooling channel is horizontal, the 
momentum conservation equation can be expressed as 
equation (9), which contains expression |ulu instead of u? 
to reflect the direction of fluid flow 
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where fis the friction coefficient and can be expressed 
as equation (10) 
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Modeling of Chemical Reaction 

In this paper, n-Decane is chosen as a fuel compound 
for the present study because it is a typical pure liquid 
hydrocarbon fuel with 10 carbon atoms in its molecule. It 


is the main component of hydrocarbon fuel and has car- 
bon number and properties similar to those of liquid hy- 
drocarbon material commonly used in aerospace applica- 
tions. 

The chemical formula of n-Decane is C,9H. and the 
critical temperature and critical pressure of n-Decane are 
617.7 K and 2.11 MPa respectively. PPD model, which 
was firstly brought out by Ward et al. in Ref. [16], is 
adopted for the thermal cracking reaction in this paper. 
The details of product distribution got from the experi- 
mental table in Ref. [16] are shown in Table 1 and the 
values in the table are average results over many experi- 
ments. The products and their distribution in a general 
PPD model are kept unchanged while the pressure and 
the temperature change at a certain pressure range. 


Table 1 Distribution of pyrolysis products of n-Decane [16] 


No. Name Formula Mass fraction 
1 Hydrogen H2 0.0021 
2 Methane CH4 0.0162 
3 Ethane C2H6 0.0267 
4 Ethylene C2H4 0.0506 
5 Propane C3H8 0.0559 
6 Propylene C3H6 0.0682 
7 n-butane C4H10 0.0419 
8 Butane C4H8 0.0777 
9 1-pentene C5H10 0.0848 
10 n-pentane C5H12 0.0630 
11 1-hexene C6H12 0.1160 
12 n-hexane C6H14 0.0538 
13 1-heptene C7H14 0.1171 
14 n-heptane C7H16 0.0506 
15 1-octene C8H16 0.1199 
16 n-octane C8H18 0.0099 
17 1-nonene C9H18 0.0468 
18 n-nonene C9H20 0.0025 


According to the general PPD model based on expe- 
rimental data used in this paper, the cracking process can 
be approximated as an overall reaction which can be ex- 
pressed as 


CoH» > 0.151, +0.143CH, + 0.256C,H4 
+0.126C,H, +0.230C,H, + 0.180C,H, 
+0. 196C, H +0.102C Hip + 0.17ICSHio (1 
+0.124C,H,, +0.195C,H,, + 0.089C;H 4 
+0.169C,H,4 + 0.072C, Hy, + 0.152C,H¢ 
+0.012C Hyg + 0.053C,H g + 0.003Cy Ho 


The cracking rate constant is generally a function of 
fuel temperature. In addition, according to Arrhenius 
expression, reaction rate k can be expressed as 


oe 
RT 
k=A-e (12) 
where A and Ea are given in Table 2 
Table 2 Reaction rate constants 
Parent fuel n-Decane 
Activation energy 
63 
Ward et al. [16] Ea, kcal/mol 
P=3.45MPa and 
T=773-873K. A factor 
2.1x1015 
A, s-l 
Activation energy 
63 
Ward et al. [9] Ea, kcal/mol 
P=3.45—11.38MPa and 
T=823-873K. A factor 
1.6x1015 
A, s-l 
Activation energy 
6442.4 


Stewart et al. [17] 
P=2.96MPa and 
T=713-808K. 


Ea, kcal/mol 


A factor 
1.0x1015.9+41.5 


A, s-l 


The fuel in the cooling channel of a scramjet engine 
works under supercritical pressure with a wide tempera- 
ture variation range, and the chemical reaction happens 
when the fuel flow through the cooling channel, increas- 
ing the complexity of fuel thermophysical properties pre- 
diction. 


Calculating method of mixture property 

Peng-Robinson equation of state (Peng & Robinson, 
1976), a Cubic Equation of State(CEOS), is used in this 
paper to predict the density of fluid in a wide range with 
a good prediction of the fuel properties near the critical 
points[18]. 

p= RT a (13) 
V-b V+2bV-b° 

where P,7,V are the pressure, the temperature and the 
specific volume. Coefficients a and b can be calculated 
using the following equations. 


a =ay{l+n0-(T/T,)°°)P (14) 
0.457247 R°T? 
ay = B (15) 
n = 0.37464 + 1.54226 — 0.269920° (16) 
RT 
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where T. and P, are the critical temperature and criti- 
cal pressure, R is the universal gas constant and œ is the 
acentric factor. The density is calculated by solving 
Eq.13, which can be regarded as a cubic equation with 
specific volume V as a variable. 

For the thermal properties of fuel, for instance, specif- 
ic heat Cp, they are calculated using relevant ideal gas 
properties and the departure functions. The departure 
function of any conceptual property F is defined as 

Pip = F geal =F (18) 


where Fidea is the value of the property as computed 
from the ideal gas relations. The departure function can 
be derived from the basic thermodynamic relations and 
the equation of state. Ideal specific heat is needed for the 
calculation of specific heat Cp. An ideal specific heat 
Cyideat is generally regarded as the function of tempera- 
ture, which can be expressed as 


C pidea /R = aT? +aT +a,T? +a;T*+a,T* (19) 


As for the viscosity and thermal conductivity of the 
fluid, the method proposed by Chung et al. is adopted"), 
And compared with the conventional method, Chung's 
method takes the high pressure and high density of the 
fluid into consideration to avoid a large relative error. 

Parameters such as M, Vc, Tc, Pc, are the molecular 
weight, critical specific volume, critical temperature, 
critical pressure and the acentric factor used for the cal- 
culation of fuel properties respectively and the parame- 
ters for the species in Eq.11 can be found in open litera- 
ture. 

The property-evaluation methods briefly described in 
this section have been used to calculate the thermody- 
namic and transport properties of n-Decane, including its 
density, specific heat, viscosity, and thermal conductivity. 
The methods have been validated using the NIST da- 
ta’! in Ref. [21]. 

As shown above, the ways to calculate the real proper- 
ties of n-Decane and its crack products have been calcu- 
lated and verified. However, the fuel will become a mix- 
ture of itself and its crack products when the crack reac- 
tion happens. So, it is very important to calculate the real 
properties of the mixture. 

A pseudo critical method is used in this paper for the 
computation of properties in a real-gas mixture and it can 


be found in Ref. [22]. According to this method, the be- 
havior and properties of a real gas mixture will be the 
same as that of a pure component, to which appropriate 
critical constants are assigned. These mixture critical con- 
stants are the functions of the mixture composition and 
pure component critical properties, and are sometimes 
called pseudo critical constants, because their values are 
generally expected to be different from the true mixture 
critical constants that may be determined experimentally. 
However, for computational purposes, they are the ap- 
propriate critical constant values for the mixture. 

According to the pseudo critical method, Eq. 13 is also 
applied for mixtures, where the critical temperature, crit- 
ical pressure, critical specific volume, and acentric factor, 
are replaced by the corresponding mixture critical con- 
stants, critical temperature, critical pressure, critical spe- 
cific volume, and acentric factor of the mixture. 

One-fluid van der Waals mixing rules”! are used for 
the mixture. According to this approach, in order to apply 
the equation of state models to the mixtures, coefficients 
a and b in Eq. 13 are replaced by the following composi- 
tion-dependent expressions: 


a => aa) (20) 
bn = >; xb; (21) 


the critical parameters of mixture can be expressed as 
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The following are the notation for Eqs. 

T n= mixture pseudo critical temperature. 

Pem= mixture pseudo critical pressure. 

Vem mixture pseudo critical molar volume. 

T.= critical temperature for component i. 

P= critical pressure for component i. 

V= critical molar volume for component i. 

x= mole fraction for component i. 

The numerical model and property-evaluation me- 
thods have been implemented into a commercial CFD 
package, Fluent, through its user coding capability. 


Numerical calculation method of equations 
All the above PDEs are solved by finite difference 


method. Firstly turned to ODEs by discretizing space 
with forward difference format, then the set of ODEs can 
be directly solved based on Simulink with suitable boun- 
dary conditions, and initial conditions. The details for 
solving the cooling channel model are shown below. 


x=L 


Fig. 2 Schematic diagram of space dimension discretization 
with boundary conditions 


Take mass conservation equation for example, the 
space dimension is divided into n equipartitions, and ge- 
nerates n+1 nodes. Point 0 is considered as a boundary 
point, and point i represents one of the middle mesh 
points. By applying a first-order difference scheme 
(backward difference scheme) into spatial derivative only, 
equation (1) can be converted to equation (13) which is a 
set of ODEs belonging to a series of mesh points 
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Similarly, by applying backward difference scheme 
into the spatial derivative of the energy conservation eq- 
uation, equation (2) can also be converted into a set of 
ODEs as equation (14) 
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Finally, applying backward difference scheme and 
forward difference scheme into term al puļu)/ ôx and 
term ép/éx of the momentum conservation equation, 


respectively, equation (9) can be converted into a set of 
ODEs as equation (15) 
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Validation of the model 


As shown in Fig. 3, an experimental set-up is built to 
investigate the flow characteristic, heat transfer and 
cracking characteristics of hydrocarbon fuel under super- 
critical conditions. 
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Fig.3 Single tube test apparatus schematic with preheater 


The facilities consist primarily of a fuel reservoir, a 


fuel pump, a fuel reactor (tube), an electric heating power 


and a fuel cooler. The tube section can easily heat fuel to 
1000 K. The test section is a GH3128 superalloy pipe. 

The pump is a constant-flux pump which can provide 
constant and continuous flow flux. The flow rate is not 
affected by the load pressure with a safety pressure as 
high as 40 MPa. The range of the flow is 0-80 ml/min. 

The test section is placing at the experimental table 
horizontally. The two ends of the pipe are wrapped with 
two copper electrodes. The pipe is heated by a direct 
current supply the power of which can be 3 kW. The 
power used to heat the pipe can be adjusted by changing 
the voltage and electric current. 

The pipe wall temperature distribution is measured 
using K-type thermocouples spot-welded to the outside 
of the pipes. In order to increase the accuracy of the 
welded thermocouple, each thermocouple was checked 
by heating the test section. The temperatures along with 
the tube are equal theoretically when the tube firstly 
heated without fuel to obtain the heat loss. The fuel tem- 
peratures at the center point of the pipe inlet and outlet 
are measured using a thermocouple inserted into the flow 
of fuel. Fuel pressure is measured with pressure gauges 
installed at the outlet of the pump and upstream back 
pressure valve, and the pressure drop is measured with 
Rosemount 3051 transducer produced by Emerson 
Company. Fuel mass flow rate is measured by MicroMo- 
tion CMF 010 mass flow meter produced by Emerson 
> Company with uncertainty less than +0.1%. The uncer- 
ge tainty associated with the measurement of wall tempera- 

ture and the temperature of fuel is estimated to be less 

z than +3 K, whereas that of pressure measurement is less 
e than 1% and that of the pressure drop measurement is 
= less than +0.075% [24, 25]. The error of the pressure 
drop is much lower than that of pressure measurement it 
is because that the sensor used to measure the pressure 
drop has a better accuracy than that used for pressure 
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(a) Inlet pressure response 


measurement (Pressure sensor produced by Micro Sensor 
Co., Ltd.) All the experimental data are recorded via a 
data acquisition system for analysis. 

The liquid and gaseous components of the fuel are se- 
parately metered and analyzed. The gas products are analyz- 
ed using a gas chromatograph/mass spectrometry (GC7900/ 
MS), which consists of flame ionization and thermal 
conductivity detectors, and the liquid portion of the 
stressed fuel is analyzed by liquid chromatography (LC). 

Two groups of experiments were conducted. The op- 
erating condition of the experiments used to validate the 
model is as follows: The back pressure is 3 MPa, the 
mass flow rate of n-Decane is 0.7 g/s, the inlet fuel tem- 
perature is 293 K and outlet fuel temperature is 773K. 
The length of the pipe is 1 m with inner diameter being 
1 mm and outer diameter being 3 mm. 

The first group of experiment is under the condition of 
inlet flow rate step change, it change from 0.4g/s to 
1g/s(Figure 4). The second group of experiment is under 
the condition of outlet throttle valve opening step change, 
it change from 50% to 55% (Figure 5). Comparison re- 
sults show the calculated data is similar to the measured 
data, the reason of dynamic characteristics will be ex- 
plained below. 


Results and discussion 


To study the effect of transient dynamic characteristics 
of cooling channel on rapidity and safety of the engine, 
two kinds of perturbation are used in the simulation. The 
input signal of mass flow rate is used to simulate the en- 
gine fuel supply process, and the change of outlet tube 
throttle coefficient is used to simulate the switching of 
flow rate valve at the outlet of the cooling channel. Dur- 
ing the calculation, the structure parameters and initial 
stable conditions are listed in Table 3. 
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Fig. 4 Comparison of calculated and measured data when inlet flow step change 
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Fig.5 Comparison of calculated and measured data when outlet throttle valve opening step change 


Table 3 The structure parameters and initial stable conditions 
of calculations 


Structure parameters 


Tube length lm 
tube wall outer diameter 4 mm 
inner diameter 2 mm 


metal density 8000 kg/m? 


metal specific heat capacity 400 J/(kg-K) 


Initial stable conditions 


outlet pressure 0.1 MPa 
heat flux 1.1 MW/m? 
Inlet mass flow rate 3 g/s 
outlet fuel temperature 750°C 


Inlet flow rate perturbation 


The inlet flow rate perturbation is from 3g/s to 3.3g/s, 
and the outlet flow rate and inlet pressure responses of 
various tube lengths, fuel temperatures and pressures are 
shown in Figure 6 to Fig. . It can be seen from that the 
longer the tube is, the slower the outlet flow rate reaches 
equilibrium state. The outlet flow rate takes almost 4 
second to reach equilibrium state in 5 meters long tube, 
compared to the 10°s of combustion, it can be concluded 
that the dynamic characteristics of engine trust control 
primarily depends on the fuel supply system. The inlet 
pressures response is shown in Figure 6 (b), inlet pres- 
sure time response curve can be divided into two parts, it 
firstly rush to a high point then slowly reduce to a lower 
stable state. it can be explained by Eq. 9 that inlet pres- 
sure increases as mass flow rate increases and inlet pres- 
sure drops as fuel velocity drops. When inlet mass flow 
rate step increases, the fuel temperature decreases slowly, 
and the fuel density drops slowly at the same time. Thus, 
the fuel velocity drops and leads to inlet pressure drops. 


It can be seen from Figure 6 (b), as the tube length in- 
crease, the inlet pressure increases more, and the over- 
shot of inlet pressure is more severe. 

As shown in Figure 7 (a), it takes the outlet flow rate 
longer time to reach equilibrium state as the outlet fuel 
temperature rises. When the outlet fuel temperature is 
300°C, the tube is full of the low compressibility liquid. 
When the fuel heated to 500°C, part of the fuel in the 
tube turns to higher compressibility super critical fluid. 
As fuel reaches 700°C, thermal cracking produces high 
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(b) Inlet pressure response 


Fig. 6 The outlet flow and inlet pressure response at different 
tube length when inlet flow step change 
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compressibility pyrolysis products, thus the outlet flow 
costs 2s to rebalance. It can be seen from Figure 7 (b), 
with the same amplitude of mass flow rate step changes, 
the higher the fuel temperature is, the slower the fuel 
temperature rebalances and thus the larger the inlet pres- 
sure increases. 

It can be seen from Figure 8 (a) that the increasing of 
fuel pressure can increase the outlet flow response rate 
feebly, because the fuel compressibility is smaller at 
higher pressure. Figure 8 (b) also shows that the overshot 
of inlet pressure is more sever at lower pressure. 
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(a) Outlet flow response 


(b) Inlet pressure response 


Fig. 7 The outlet flow and inlet pressure response at different 
outlet fuel temperature when inlet flow step change 


Outlet throttle valve opening perturbation 


As the outlet throttle valve opening decreases by step, 
the fuel in the tube rush out of outlet then the flow rate 
recovers to the initial state because the inlet flow rate 
keeps unchanged. In section 4.1, it is known that when 
inlet flow changed, the outlet flow response time is very 
long. It can be seen from Figure 9 (a) that controlling the 
outlet valve opening can be a method to increase the re- 
sponse speed of the engine fuel supply system. Figure 9 
(c) and Figure 9 (d) shows the change of the outlet throt- 
tle valve opening has a few effects on the wall tempera- 
ture and fuel temperature. 
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Fig. 8 The outlet flow and inlet pressure response at different 
pressure when inlet flow step change 


Conclusion 


To investigate the transient dynamic characteristics of 
cooling channel during the process of fuel flow supply 
and the change of cooling channel outlet valve opening, 
this article focuses on the simulation of a one-dimen- 
sional model with high temperature hydrocarbon fuel at 
pyrolysis condition. Tube length, working pressure and 
heat load are chosen to investigate the effect on channel 
parameters dynamic characteristics. 

An important result was that the time constant of out- 
let flow is 2 to 5 seconds, and time constant increased 
with increasing tube length and heat load, working pres- 
sure has few effects on it. As inlet fuel step changed, the 
inlet pressure appeared overshoot, thus caused transient 
impulse on pumping power. Increasing the tube length, 
heat load and decreasing the working pressure all lead to 
higher transient impulse. To increase the response speed 
of the engine fuel supply system, the step decrease of 
outlet throttle valve opening was a good way without 
large fluctuation of fuel and channel wall temperature. 
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(d) Outlet fuel temperature response 


Fig. 9 The outlet flow, temperatures and inlet pressure response at different valve change amplitude when outlet throttle valve 


opening change by step. 
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